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The transverse broadening of jets produced in deep- inelastic scattering (DIS) ofT a large nucleus 
is studied in the coUinear limit. A class of medium enhanced higher twist corrections are re-summed 
to calculate the transverse momentum distribution of the produced coUinear jet. In contrast to 
previous approaches, re-summation of the leading length enhanced higher twist corrections is shown 
to lead to a two dimensional diffusion equation for the transverse momentum of the propagating 
jet. Results for the average transverse momentum obtained from this approach are then compared 
to the broadening expected from a classical Langevin analysis for the propagation of the jet under 
the action of the fluctuating color Lorentz force inside the nucleons. The set of approximations 
that lead to identical results from the two approaches are outlined. The relationship between the 
momentum diffusion constant D and the transport coefficient q is explicitly derived. 



PACS numbers: 12.38.Cy, 25.75.-q, 13.60.Hb, 13.87.-a 



I. INTRODUCTION 

The study of dense matter through its effect on hard 
partonic jets is now an established science. Experiments 
at the Relativistic Heavy-Ion collider (RHIC) and the 
upcoming Large Hadron Collider (LHC) have consider- 
able resources dedicated to the study of the modification 
of high transverse momentum particles produced in the 
fragmentation of partonic jets which originated within 
and propagated through the produced dense matter 
There currently exists a large number of theoretical mod- 
els which provide quantitative estimates of the observed 
data 

While the different models have rather diverse origins, 
within the approximations made, similarities are perva- 
sive. A well know example is the quantitative similarity 
between the schemes of Refs. Q and as demonstrated 
in Ref. The object of this article is two fold: primar- 
ily, the aim is a re-derivation of the transverse broad- 
ening of a hard jet brought about by its propagation 
in dense matter with fluctuating color fields. Secondly, 
this represents the first effort to find a physical resem- 
blance between the Higher- Twist and Finite temperature 
field theory jet-quenching schemes Refs. [1, The fo- 
cus will lie predominantly on the Higher-Twist approach 
of Ref. [J. In this scheme of energy loss, one tends to 
re-sum a class of diagrams which encapsulate the lead- 
ing medium length enhanced power corrections to the 
basic process of a hard parton radiating a soft coUinear 
gluon. Such gluon showers, emanating from hard partons 
in vacuum, often referred to as vacuum energy loss ^ (or 
zeroth-order energy loss lead to the DGLAP evo- 
lution of parton-hadron fragmentation and structure 
functions. Such showers, in a medium, are influenced 
by the time (or length) dependent transverse momentum 
fluctuations experienced by the hard propagating jet and 
its radiated gluons and lead to a loss of the energy of the 
propagating parton. The effect of the transverse "kicks" 
from the medium on the systematics of energy loss may 
be divided into two parts. The effect on the propagation 



of the jet as well as its radiated gluons and the effect on 
the radiation vertex. The current work, intends to focus 
exclusively on the former i.e., the effect of the color fields 
of the medium on the transverse momentum fluctuations 
or broadening experienced by a single hard parton as it 
traverses the medium. Radiation and energy loss will be 
dealt with in a future effort. 

It is well known that hard partons, traversing dense 
matter, tend to pick up a transverse momentum which 
depends on the length of the medium The calcula- 
tion of this length dependent broadening is achieved via 
the re-summation of a class of higher twist diagrams [ll| 
which are enhanced by the length traversed by the jet 
in the medium. Within the assumptions made in such a 
calculation, the broadening may also be obtained from a 
purely classical analysis of a hard charged particle mov- 
ing under the influence of a color Lorentz force. In this 
sense, there exists a physical similarity with the basic 
picture underlying the jet quenching scheme of Ref. Q 
which also admits a kinetic theory description of hard 
partons moving under the influence of soft fields. 

However, the existing literature on the inclusion of 
power corrections to energy loss or jet broadening pro- 
cesses has yet to yield such a simple physical picture. 
This will be the object of the current article. This is 
hardly the first attempt to identify and re-sum the class 
of higher twist corrections required for the computation 
of jet broadening in extended media. A thorough deriva- 
tion is provided to elucidate the basic assumptions and 
approximations made. Similar to previous derivations, 
the original parton emanating from the hard scattering 
is assumed to posses a vanishing transverse momentum 
distribution ((5^(p_|^)). The derivation is carried out in the 
high energy limit; as a result, the coupling of the hard 
parton with the medium is assumed to be weak. A deriva- 
tion of the transverse broadening, when the parton cou- 
ples strongly with the medium (albeit for a heavy quark) 
has been presented in Ref. p^ . In contrast to previous 
derivations, re-summation of the leading higher twist cor- 
rections, leads to a diffusion equation for the transverse 



momentum distribution of the final parton on exit from 
the dense medium. The diffusion equation is solved and 
the relationship between the diffusion tensor and the en- 
ergy loss transport coefficient q derived. The compu- 
tation of the first non-vanishing moment of transverse 
momentum distribution {p^) is then re-derived within a 
classical approximation where a colored charge moves un- 
der the influence of the fluctuating Lorentz force within 
nucleons. This calculation is carried out using a simple 
Langevin analysis. The similarities between the two re- 
sults and its implications are discussed. 

The paper is organized as follows: in the next sec- 
tion, the leading twist parton distributions will be fac- 
torized at leading order. In Sect. Ill, the class of higher 
twist diagrams which arc length enhanced at a given or- 
der m will be identified and the leading contributions 
to the hadronic tensor will be calculated. In Sec. IV, 
the hadronic tensor will be factorized into hard and soft 
piece and the m"^ transverse momentum derivative of 
the hard part computed. In Sec. V, an all order re- 
summation is carried out and the transverse momentum 
diffusion equation is derived. The diffusion equation will 
be solved and the moments of the transverse momentum 
distribution computed and related with the energy loss 
parameter q. In Sect. VI, the Langevin analysis is carried 
out. Concluding discussions and future directions will be 
presented in Sect. VII. 



II. LEADING TWIST AND PARTON 
DISTRIBUTION FUNCTIONS 

The focus of this article is restricted to the semi- 
inclusive process of DIS off a nucleus in the Breit frame 
where one jet with a transverse momentum l± is pro- 
duced, 



e(Li) + A{p) e{L2) + J(/±) + X. 



(1) 



In the above equation, Li and L2 represent the momen- 
tum of the incoming and outgoing leptons. The incom- 
ing nucleus of atomic mass A is endowed with a mo- 
mentum Ap. In the final state, all hadrons [hi, h2, ■■■) 
with momenta pi,p2, ■ ■ ■ are detected and their momen- 
tum summed to obtain the jet momentum and X denotes 
that the process is semi-inclusive. 

The kinematics is defined in the Breit frame as 
sketched in Fig. [1] In such a frame, the virtual photon 
7* and the nucleus have momentum four vectors q. Pa 
given as, 



q = L2-L1 = 



2q- 



,g-,0,0 



Pa = A[p+, 0,0,0]. 



In this frame, the Bjorken variable is obtained as xb — 
QV2p+q-. 





FIG. 1: The Lorentz frame of the process where a nucleon in 
a large nucleus is struck by a hard space-like photon. 



The differential cross section of the semi inclusive pro- 
cess with a jet with transverse momentum l±_ and four- 
momentum I may be expressed as 



(PL2(P-h 



EM T 



2tisQ 



~d?T7' 



(2) 



where s = [p + ii)^ is the total invariant mass of the 
lepton nucleon system. The reader may have already 
surmised the form of the leptonic tensor as, 



(3) 



In the notation used in this paper, \A;p) represents the 
initial state of an incoming nucleus with A nucleons with 
a momentum p per nucleon. The general final hadronic 
or partonic state is defined as \X). As a result, the semi- 
inclusive hadronic tensor may be defined as 



W^'''=^2^^)b\q + PA-Px) 



X 



X {A;p\J^mX){X\r{0)\A;p) 



= 2Im 



d'ye^'>-y{A;p\J^{v)rmA;p) 



(4) 



where the sum (X^jc) runs over all possible hadronic 
states and is the hadronic electromagnetic current 
(J^ ~ Qq'tpq^'^^q, whcrc Qq is the charge of a quark of 
flavor q in units of the positron charge e). It is under- 
stood that the quark operators are written in the inter- 
action picture, and factors of the electromagnetic cou- 
pling constant have already been extracted and included 
in Eq. ([2]) . The leptonic tensor will not be discussed fur- 
ther. The focus in the remaining shall lie exclusively on 
the hadronic tensor. This tensor will be expanded order 
by order in a partonic basis and leading twist and maxi- 
mally length enhanced higher twist contributions will be 
isolated. 

The leading twist contribution is obtained by expand- 
ing the products of currents at leading order. This contri- 
bution may be expressed diagrammatically as Fig. [2] and 
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quantitatively expressed as (we also take the average over 
initial states and sum over final states to obtain) 



a 



2'KX} 



(5) 



In the above equation, expresses the probability to 
find a nucleon state inside a nucleus with A nucleons and 
fqixs) is the quark structure function within a single 
nucleon state, this represents the expectation value of a 
twist 2 operator: the quark number operator. Through- 
out this paper, the light-cone component notation for 
four vectors [p = [p'^,p~ ,p±\) will be used, where. 



P 



P + Pz 



p =p 



Pz 



(6) 



In the coUinear limit, the incoming parton is assumed to 
be endowed with very high forward momentum (po^ 
xqp^.Pq 0) with negligible transverse momentum 
Po± *S Po ■ Within the kinematics chosen, the pho- 
ton also has no transverse momentum. As a result, the 
produced final state parton also has a vanishingly small 
transverse momentum [i.e., with a distribution ^^(p^)). 
As a result, the transverse momentum distribution of the 
produced parton is obtained as the differential hadronic 
tensor, 



(7) 



as the cut line in Fig. [2] would indicate. The outgoing 
quark will eventually fragment into a jet of hadrons and 
thus has a positive virtuality which is dependent on 
the scale of the hard scattering. In a high energy 
scattering event Aqcd << m < Q, which allows futher 
re-scatterings of the quark to be treated perturbatively, 
yet m « q~ which leads to its identification as a single 
jet. The scale m may also be used to define the scale of 
the final fragmentation function. In the remaining sec- 
tions, both the virtuality of the produced quark and its 
initial transverse momentum distribution will be taken 
to be vanishingly small. These approximations are car- 
ried out in the interest of simplicity and both quantities 
should be understood to be present. 



III. HIGHER TWIST AND TRANSVERSE 
BROADENING 



Higher twist contributions are obtained from diagrams 
which include expectation values of more partonic opera- 
tors in the medium [isj e.^., the gluon field strength oper- 
ator product F^'^ {y)F^ (Q). To obtain higher twist con- 
tributions, higher orders need to be included. A diagram 
with 2n gluon insertions may contribute to twist m < 2n. 
Issues relating to the generalized factorization [3, ES] of 
such contributions will not be dealt with in this effort; 
the focus will be to obtain an effective description of the 
propagation of a hard parton in transverse color fields. 
In all calculations, the high energy and hence small g 
limit will be assumed, as a result all diagrams containing 
the four-gluon vertex as in the left panel of Fig. \3\ are 
suppressed compared to diagrams with the same number 
of gluon insertions on the hard line which are directly 
connected to the soft matrix elements. As a result, all 
such diagrams will be ignored. 





FIG. 2: Tlie Lowest order and leading twist contribution to 

It should be pointed out that in DIS on a nucleon, the 
final transverse momentum distribution of the produced 
jet is never a strict 5-function but rather a Gaussian with 
a width that depends on the scale and xb of the scatter- 
ing. Neither is the final outgoing quark strictly on shell 



FIG. 3: A higlier order contribution suppressed by due to 
four gluon coupling. 

The aim is to isolate the higher twist contributions at 
a given order and twist that carry the largest multiple 
of length L A^^^. In the language of power count- 
ing, one looks at the combination a™L" (usually n < m) 
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and focuses on diagrams with the maximum n. The mul- 
tiples of L are obtained by insisting on conditions that 
lead to the largest number of propagators going close to 
their on-shell conditions. This is easily achieved by the 
hard parton having the maximum number of space like 
exchanges with the medium. Such diagrams may or may 
not include radiated gluons. The introduction of radiated 
gluons leads to the calculation of energy loss of the parent 
parton Computation of diagrams without radiated 
gluons deal with the propagation of the parent parton 
without radiative energy loss. Diagrams where a single 
gluon line originates in the medium and then splits into 
two prior to attaching with the hard parton (left panel 
of Fig. |4]) represent virtual corrections to radiative dia- 
grams and fulfill unitarity conservation at higher order. 
While radiative contributions usually require at least one 
propagator to be off-shell, they are enhanced by large log- 
arithms and thus represent a somewhat different power 
counting. Such diagrams along with their real counter- 
parts will be dealt with in a future publication.. 

Diagrams where multiple gluon lines fuse into a sin- 
gle gluon prior to attaching with the hard parton (right 
panel of Fig.d]) do not produce the same length enhance- 
ment as that of the individual gluons attaching directly 
to the hard parton. The reason behind this is that if the 
intermediate gluon line is on-shell then it forces the next 
quark propagator to go off shell or vice-versa. In cither 
case, it is suppressed by one factor of L (this is treated 
in more detail in Appendix A) . There is an exception to 
this condition that occurs in the region where the forward 
momentum fraction of the gluon x is very small. In such 
cases, due to the saturation mechanism, soft gluon pop- 
ulations may be enhanced [3| ^nd such diagrams may 
indeed produce considerable contributions [ij, [11] . The 
current analysis will be assumed to be carried out out- 
side the saturation regime. As a result, such gluon fusion 
contributions will be ignored. 



Note, there is no particular ordering of the gluon lines. 
Also, the gluon lines are not meant as propagators but 
rather as field insertions at a point y^; hence there is no 
meaning associated with crossed gluon lines. The en- 
tire set of n -|- n' vertex insertions (with the gluon vec- 
tor potentials contracted with the nucleus) may then be 
connected by quark propagators in {n + n')! ways (this 
overall combinatoric factor is removed by the {n + n')! 
which appears in the denominator from the perturba- 
tion expansion). The reader may question the focus on 
this sub-class of possible diagrams at order n + n', where 
all gluon field operators are contracted with the nuclear 
state. Diagrams of the same order, where the gluon field 
operators are contracted with each other will contain cut 
(or uncut) gluon lines, and these will represent real (or 
virtual) contributions to gluon radiation diagrams. A sub 
class of contributions such as those in Fig. |4] may indeed 
be included by a redefinition of the effective gluon vector 
potential. As such, they would represent small additive 
contributions which do not influence the leading length 
enhanced behavior of the diagrams of Fig. [5] and will not 
be discussed further. 





FIG. 5: An order 2?! contribution to W^". This contributes 
to twist m < 2n. 



FIG. 4: Left panel shows a virtual contribution to a radiative 
correction. Right panel is a gluon fusion contribution, which 
is suppressed outside the saturation region. 

The diagrams under discussion have the general form 
of Fig. [5l The ellipses in between gluon lines in Fig. O 
are meant to indicate an arbitrary number of insertions. 



At the [n + n') order, we envisage a general contribu- 
tion where there are n incoming gluon lines originating in 
the locations yi, . . . ,?/„ in the amplitude, whereas there 
are n' gluon lines at locations yj^ , . . . , y^, in the complex 
conjugate. The cut line is generically labelled as I. The 
Feynman rule for such a contribution to the hadronic 
tensor is given as, 
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qf - ie 



X e 



(8) 



The momenta, 'i'ii'/j label the momenta of the various 
fermion lines. Note, there is overall momentum eonser- 
vation within the diagram. As a result, a shift in all 
momenta may be executed as a consequence of momen- 
tum conservation at each vertex: qi ^ X]m=o P™ ^ 9 ^^"^ 
q'j = J2k=^o P'k + 1' where Pi,p'j are the momenta brought 
in by the i^^ line in the amplitude and the j"^ line in 
the complex conjugate. It should be pointed out that all 
such Pi,p'j integrals run from — cx) to oo and thus may 
indicate both incoming and outgoing momenta. Using 
these substitutions, the phase factors in Eq. ^ may be 
written as 



n— 1 n —1 

X exp -jy„ • |/ - + 



X exp 



n'-l 



(9) 



3=0 



where, it is understood that ?/q is the origin and = y. 
At this point an n}^ momentum may be introduced, via 



1 = J dSnS' [l p^Pk q^ 



(10) 



This leads to a considerable simplification of the phase 
factor as 



r — exp 



n n — 1 

i=0 j=0 



— 1 

■0 j=0 



(11) 



Note the complete absence of the hard photon momen- 
tum q from the phase factor. 



The approximations stemming from coUinear dynam- 
ics may now be instituted. The calculation is carried out 
in the Breit frame at very high energy. As a result, all 
momentum lines that originate in the target are domi- 
nated by the large -I- components of their momentum, 
followed by their transverse coordinates, i.e., 



Pi »Pt± »Pr 



(12) 



In most cases the above condition will allow us to prac- 
tically drop all — components of momentum from the 
expression for the hadronic tensor and focus solely on 
the -|- and _L components. This procedure will be carried 
out not only for the momenta but also for the field op- 
erators. In this effort, calculations will be carried out in 
the covariant gauge. In covariant gauge calculations in 
the Breit frame at very high energy, the dominant com- 
ponents of the vector potential are the forward or (-I-) 
components [3, i.e., A'' ~ g-''A+. As there are no 
cancellations involving the vector fields, there will arise 
no further need to consider the _L components. In the 
corresponding calculation in light-cone gauge ~ 
and the largest components are the A-^ components. It 
should be pointed out that the — components arc only 
being dropped from locations where they appear in ad- 
dition to the larger -l-,-L components, i.e., they are not 
dropped from the phase factors. The coUinear approxi- 
mation may also be used to simplify the quark operators, 
as 



+ i^my)i±4m)] 



(13) 



where, once again due the high energy limit, all terms 
coUinear with the + direction are enhanced by boost com- 
pared to the _L direction and thus the second term above 
maybe dropped. 

The hadronic tensor of Eq. ([5]) may be expressed as 
the formal convolution of four terms. 



d2;^Dy2)pT(p)r(p,y)A/(y) 



(14) 



where, T{p) denotes the pure momentum component 
of the integrand, Af (y) is the pure position dependent 
multi-operator matrix element and r(p,y) is the phase 
factor of Eq. (jlip that convolutes positions and momenta. 
The bold face quantities p,y represent an array of mo- 
menta and positions. 



p = bo, ■ • ■ ,Pn;p'o, ■ ■ ■ ,p'n'-i;l], 



6 



The integration measures and Dy denote a product 
of integrals over the different four vectors contained in 
the arrays above. In Eq. ((M]), the integrations over the 
Ughtcone components of the cut line I have been per- 
formed using two of the four delta functions introduced 
in Eq. 1101 The remaining two components of the trans- 
verse momentum integration are 



dH^5'^{l^^- K^_) = 




M{y) = TrlA-p 



A-p 



(17) 



While, the phase factor which convolutcs both these 
terms is essentially contained in Eq. pT|) . 

The complete absence of the (— ) components of the 
momentum, from all expressions except for the phase fac- 
tors allows for the p~ and p'~ integrations to be done, 
resulting in the localization of the process on the negative 
light-cone, i.e., 



where, K±^ is a representative of the sum of the trans- 
verse momenta brought in by the n gluon insertions. The 
solely momentum dependent terms in the integrand, post 
coUinear approximation, may be expressed as, 



T{p) = Tr 



7 7'^ 



n (po+g^) + 7 g -7 -Po 



X 7 



X 7 



7' ( T.Pt + <l' 

.1=0 



7 ' g - 7^ ■ ( E W 

,i=0 



/n-1 

2 E P+ ) <?- - 0' 

\i=0 



+ +i^q -7 



2T:5l2q- [Y^pt + q 




X 7 



,'i=0 



7- ( e'p'; + 9+ ) + ^^9- - 7^ • [Ylp'j 



7 (p'o+g )+7 g -7 -P'o u 



i=0 j"=0 " 

Z; ip'Jiy'j-y'^,) 

X e 



J=0 



i=0 



In the above equation, ?/'q = and as may be noted from 
the definition of the location arrays is not being inte- 
grated over. As a result, this constrains all the negative 
light cone locations in the above equation to the origin. 

The terms in the numerators of Eq. p6p may be sim- 
plified with observation that 7^7^ = = 0, as a 
result only the 7'*' and 7^ terms survive the spin sum. 
In the extremely high energy limit where q~ 3> P^P'f ■ 
terms proportional to the transverse components of the T 
matrices may also be ignored. The terms in the denomi- 
nators may be simplified further with the replacements 



= 2xBP^q ; pi = Xip~^ ; p'^ = 2:'p+ 



(16) 



k=0 



(19) 
(20) 

(21) 



(=0 



which is the short distance momentum dependent part 
of the integrand. The ellipses in the above equation are 
meant to indicate the presence of propagators with a 
growing number of additive momentum factors brought 
in by the gluon insertions. Note that all appearances 
of the (— )components of the momenta (p^) have been 
ignored in the above equation. The long distance non- 
perturbative position dependent part of the integrand is 
given as 



In the above equation, i, the index of the unprimcd mo- 
menta (both for the longitudinal and transverse compo- 
nents) runs from to n, whereas, the index of the 
primed momenta runs from to n' — 1, i.e., one less than 
maximum. 

With the above simplifications, the momentum depen- 
dent part of the hadronic tensor W'^'^ , as sketched in 
Eq. pB]) assumes the form. 
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T{p) = 4 



(2p+)" 



+ri' + l 




(22) 



The integrals over the momenta p^,p'^ may be re- 
expressed in terms of momentum fractions, i.e., dp^ = 
p'^dxi- The n + n' + 1 integrals over (+) momenta in 'Dp 

lead to an overall factor of , which is cancelled 

by the similar factor appearing in the denominator of 
Eq. ([22|l . As a result, the integral measure now has the 
appearance. 



r+ = exp 



-iixB + Y^ x'o P^iVn - y'n ) 



2 = 

n'-l 

X Y[ 6xp 



ix'jp+{y'j ~2/'„') 



(26) 



where, the second line in the equation above involves 
only momentum fractions and locations from the left- 
hand side of the cut, where as the last line involves mo- 
mentum fractions and locations from the right hand side 
of the cut line. The integrations over the remaining lon- 
gitudinal momentum fractions, may now be performed 
starting from the propagators adjacent to the cut and 
proceeding to the propagators adjacent to the photon 
vertices. 

The first such integration, involves the propagator 
from the 2nd line of Eq. (P^ . Isolating the piece that 
depends on the fraction Xn~i yields the integral, which 
may be performed by closing the contour of x^_i with 
a counterclockwise semi-circle in the positive imaginary 
direction, 



i=0 



A dx. d^p.^ -yr' dx'^ rf^p;.^ 



27r (27r)2 J-l 27r (27r)2 ' 

i=0 ^ ' j=Q ^ ' 



The remnant phase factor may be decomposed into a 
longitudinal piece and a transverse piece, which depends 
on the corresponding components of the momentum ap- 
pearing as arguments, i.e., 

n n —1 

n n' — 1 

j=o 



X lie"'" '•"^ ^'^'Y[e-'Pi ■^y.-y^) (24) 

i=0 



The cut line, in the diagram of Fig. [5] is indicated by 
the delta function appearing in Eq. (j22p . Integrating over 
the last momentum fraction Xn with the aid of the delta 
function leads to the condition that 



n n—1 

1=0 i=0 



(25) 



Instituting this condition leads to the separation of the 
longitudinal phase factor r+ into a left and right piece. 



dXn-l 



2tt 



n-2 



„n-l 



= iSiVn - Vn-l) 

Ti-2 
i = 



xb - le 



(27) 



The effect of performing the above integration is the 
incorporation of both the real and imaginary parts of 
the above propagator into the overall expression of the 
hadronic tensor. It has the physical effect of propagat- 
ing the quark from y~_i to y~ . Similarly, the integration 
over the propagator to the immediate right of the cut line 
may be carried out by closing the contour of with a 

clockwise semi-circle in the negative imaginary direction: 



dr' 



■ / + / ' — ' — \ 



2n n'-2 ^ ^ 

j=0 



Xb + «e 



3=0 



+ ( I — I — \ 

p iv „i_i-y „/) 



(28) 



Incorporation of the results of the above two integrals 
into the longitudinal phase factors leads to the expres- 
sion. 
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r+ = cxp 



■ 'n — , - in -I-/ — 

-iXoP Vn +T'X dP y n 



-IXbV y.n-1 +tXBP y n'-l 
n'-l 



i=0 



J=0 



Tl-2 



X 



Yl exp [-ix,p+{y^ - 



i=0 
n'-2 

X exp 



(29) 



where, overall momentum conservation was invoked to 
define the new variable 



n'-l 



(30) 



In the above expression, x]j for < z < rt — 1 are defined 
in Eq. ([^D|) . and x'^^j for < j < n' — 1 are defined in 
Eq. (PT|) . The n"^ transverse fraction is set by the 
delta function arising from the cut line and is given as in 
Eq. 

The trend of the longitudinal momentum fraction in- 
tegrals is now clear and the general result for the integra- 
tions over the remnant momentum fractions in the phase 
factor may be carried out. The full hadronic tensor of 
Eq. ([5]), post integration over all longitudinal fractions, 
may now be expressed as 



,2, n n 

^Udy-dY,l[dy'-dYi 



<Pp\ -i-r d^p'- 



ii (27r)2 U (2^) 



i=0 



X e 



i=0 



X e'" -D'' " J e~ 
3=0 

1 1 



X 



n%r-yr-i) n -yVi) 

i=n. j—n' 



X (yl;p|^(2/-,yi)7+V;(0)Tr 



i=l 



X llt^^A+{y';,y,i 



(31) 



The expression derived above is completely general, in 
the sense that no assumption regarding the nature of the 
nuclear state has been made. In the next section, a fac- 
torization of the above hadronic tensor into a part that 
is solely dependent on hard momenta and a part depen- 
dent on soft momenta will be carried out and simplifying 
assumptions regarding the nuclear state made. 



IV. FACTORIZATION AND GRADIENT 
EXPANSION 



Up to this point, no approximation regarding the na- 
ture of the state \A;p) or the action of the quark and 
gluon operators on this state has been made. We now 
approximate the nucleus as a weakly interacting homo- 
geneous gas of nucleons. Such an approximation is only 
sensible at very high energy, where, due to time dilation, 
the nucleons appear to travel in straight lines almost in- 
dependent of each other over the interval of the interac- 
tion of the hard probe. In a sense, all forms of correlators 
between nucleons (spin, momentum, etc.) are assumed 
to be rather suppressed. As a result, the expectation 
of the n + n' + 2 operators in the nuclear state may be 
decomposed as 



n-t-n 

{A;pmy-,y^)j+ij{0) J] Al{y,)\a;p) 

n+n 

= AC^^{p,my-,y^h+m n 

1=1 

+ Cplp.(Pil'^(2/",y-i-)7+^(0)bi) 



n+n' 

X {P2\ n Kiy^)\p2) 



(32) 



where, the factor C^^ represents the probability to find 
a nucleon in the vicinity of the location y, which is a 
number of order unity (it is the probability to find one 
of A nucleons distributed in a volume of size cA within 
a nucleon size sphere centered at y). The remaining co- 
efficients Cjf^ represent the weak position correlations 
between different nucleons. The overall factor of A arises 
from the determination of the origin (the location in the 
equation above) in the nucleus, which may be situated on 
any of the A nucleons. 

It is clear from the above decomposition that the 
largest contribution arises from the term where the ex- 
pectation of each partonic operator is evaluated in sepa- 
rate nucleon states as the yi integrations may be carried 
out over the nuclear volume. As a nucleon is a color 
singlet, any combination of quark or gluon field strength 
insertions in a nucleon state must itself be restricted to 
a color singlet combination. As a result, the expectation 
of single partonic operators in nucleon states is vanish- 
ing. The first (and hence largest) non-zero contribution 
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emanates from the terms where the quark operators in 
the singlet color combination are evaluated in a nucleon 
state and the n + n' gluons are divided into pairs of sin- 
glet combinations, with each singlet pair evaluated in a 
separate nucleon state. This requires that n + n' is even 
and may lead to a maximum overall factor of 



a 



Pl,P2, 



^[(n+n'+2)/2] 



(33) 



in the large A limit. It should be pointed out that large 
contributions may also arise, in principle, when n + n' is 
odd. In this case, the two quarks and a gluon arc consid- 
ered in the singlet combination with the remaining glu- 
ons evaluated in singlet pairs in the remaining nucleons. 
Here we institute the experimental observation that (n)- 
parton observables are much smaller than (n — l)-parton 
observables. This is only true, once again, outside the 
saturation regime, as discussed at the beginning of this 
section. Such odd parton expectations also become im- 
portant in cases where nucleon polarization effects are 
being studied where the spins of the nucleons are not 
averaged over [l^. In this effort, the focus remains ex- 
clusively outside such regions, as a result we ignore all 
terms with more than two quarks or two gluons per nu- 
cleon. 

Further simplifications arise in the evaluation of gluon 
pairs in a singlet combination in the nucleon states by 
carrying out the y± integrations. The basic object under 
consideration is (ignoring the longitudinal positions and 
color indices on the vector potentials) 



d'yld'y"Ap\A+ifi_)A+iy'l)\p) 



X e^P--yHp\A+{yj2)A+{-yj2)\p), 



(34) 



where, y± is the transverse gap between the two gluon in- 
sertions and p± = (p^ +p-'j)/2. The physics of the above 
equation is essentially the transverse translation symme- 
try of the two gluon correlator in a very large nucleus. 
One will note that the two dimensional delta function 
over the transverse momenta has removed an integration 
over the transverse area of the nucleus thus reducing the 
overall A enhancement that may be obtained. This is 
then used to equate the transverse momenta emanating 
from the two gluon insertions in the amplitude and com- 
plex conjugate amplitude. This also simplifies the lon- 
gitudinal phase factors which now depends solely on the 
difference of the longitudinal positions of the two gluon 
insertions. 

Up to this point, the collinear approximation has been 
used to simplify the expressions for the hadronic tensor, 
without the introduction of factorization. The separa- 
tion of the hadronic tensor into a hard short distance 



piece and a soft long distance contribution may now be 
accomplished. All factors in Eq. ([3T|) which contain the 



hard scales p'^,q~ constitute the hard part. All factors 
that depend solely on the soft _L momenta and distances 
along with the matrix clement constitute the long dis- 
tance element, i.e., the first part of the integrand in the 
3'*^ and 4*'' lines in Eq. ([3T|) belongs in the hard part 
along with all phase factors which contain a factor cc^p"*" 
or x'''jjp~^ as part of their arguments. The purely trans- 
verse phase factors such as e'Kp[ip± ■ y±] belong in the 
soft part along with the matrix elements. Thus, we may 
decompose the hadronic tensor as, 



Vy^p^H{p+,q-,p^,y)S{p^,y). (35) 



In the above equation, y and p± are representative of the 
entire set of distances and transverse momentum that ap- 
pear in Eq. ()31|) . One now invokes the collinear approxi- 
mation in expanding the hard part as a Taylor expansion 
in transverse momenta around the origin pY ^ 0- the 
case of a symmetric cut with n = n' there are 2n gluon 
insertions and as a result, as many derivatives, which in- 
volve n different transverse momenta. All terms of the 
form 



2 dp\"dp\ 



H 



p^=oPl PI 



where, m < n, yield gauge corrections for the contri- 
butions with 2m gluon insertions and as many deriva- 
tives [l^. The genuine 2n twist correction at this order 
has m = n. Expanding to this order, one obtains the 
generic term. 



n 



2 dp\°'dp'' 



-H{p+,q ,pYJ|p^=o 



X p^^"p\Ul{yr,yl/2)A+, {y' ^ , -yl/2). (36) 

Where, we have assumed the result of Eq. ([M]) and rein- 
troduced the color indices and longitudinal locations. 

Using integration by parts over the transverse distance 
y\^ one may convert the product AJ. (i/^", yY/2) — > 

■if?" (^i" 1 2/1/2) ■ In the extreme collinear limit, in 
the presence of a hard scale such that g is small, one 
may make the approximation. 



d^A^„ 



F„ 



(37) 



where, i^aj represents the gluon field strength. Carrying 
this out consistently on the two gluon operator in the 
nucleon state of Eq. and ignoring derivatives of the 
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field strength 



0] we obtain, 



d'y±p\_"p\e^P--y-{p\A+{y^/2)A+{-ijj2)\p) 
d^y^e'P^-y^l{p\F+^{y^/2)F+f{-y^/2)\p) 



d'y±e^P--y-^{p\F+>'iyj2)F+i~yj2)\p)m 



In the last line of the above equation we have averaged 
over the spins in the two field strength expcction in the 
nucleon state with the constraint that the operator being 
evaluated in the nucleon be a spin singlet. The nucleon 
states are always assumed to be spin singlets or in spin 
averaged states. 

The hard part which consists of 2n transverse momen- 
tum derivatives may now be simplified further 



n- - 

i dp\ 



5\li_ - K\) = {Viy 5\li_ - R\). (39) 



With the derivative expansion (at vanishing transverse 
momenta p\ — > 0) imposed on the hard part _ff , it no 
longer has any functional dependence on the transverse 
momenta. The integrations over the transverse momenta 
may now be included completely into the soft part. The 
action of the derivatives remains on the transverse 6- 
function as the only non- vanishing contribution from the 
colhnear expansion on Eq. pip . 

In the interest of simplicity, we have considered the 
symmetric case of gluon insertions above, where n = n' . 
The other cases may be easily computed in similar fash- 
ion (the case where n = n' — 2 is considered in Appendix 
B). The longitudinal integrals, due to color confinement 
yield the requirement that the longitudinal locations of 
the two gluons which act on the same nucleon state be 
in close proximity. One now tries to identify the most 
length enhanced term by isolating the maximum number 
of unconstrained dy^ integrals. Note that, due to color 
confinement (ignoring color and spin indices). 



dy dy' d^y 



r-ixnP'^iy ~v' ) (AO) 

(2^)2^ ^ ^ 



X e'P^-y^{p\Fiy- ,yj2)F{y'- ,~yj2)\p) 
= ldY-JdSy-d^^^e-^r>-Sy- 

X e^P^-y^PlFiSy /2,y^/2)F{-6y- /2,-y^/2)\p) 



where, xd 



2p+q- 



is a function of the pj^ itself. 



In the equation above, all three integrals over 5y~ ^ y± 
are limited by the nucleon size. However, the integral 
over Y~ ^ under the assumption of longitudinal transla- 
tion invariance of the two point correlator, may span the 
entire length of the nucleus. Each such integral yields a 
factor of L~ ^ A^^^ from the unconstrained integra- 
tion. The equating of the pairs of transverse momenta 



that appear in each two-gluon correlation, as well as the 
relation between the longitudinal momenta from the 9- 
functions in Eq. pip , require that the largest transverse 
momentum broadening and largest length enhancement 
arises from the terms where the gluon correlations are 
built up in a mirror symmetric fashion, i.e., where the 
gluon insertion at is contracted with that at y''. This 
reduces the focus on terms where n = n' which produce 
a transverse momentum broadening of order A"/'^ and 
provides an aposteriori justification for our discussion of 
symmetric terms. 

One may average colors of the quark and gluon field 
operators. 



srab 



{p\F^F'\p) = 



[N^ 1) 



{p\F''F-\p) 



(41) 



(p|Vi.7+V'jb> = j^{p\i^i^^\p)- (42) 



This reduces the overall trace over color factors to 



1 



NciN^ - 1)" 



Tr 



i—l j—n 



[N^ - 1)" i2NcY' 



(43) 



If the original parton was a gluon, one would simply re- 
place Cp with Ca in the above equation. 

The remaining n longitudinal position integrals for the 
gluon insertions may be simplified as 



ff[dyrO{y--y-_,) = l-ff[dy,^ 



(44) 



Invoking the above simplifications, the leading length 
enhanced contribution at order 2n to the differential 
hadronic tensor is obtained as (in the following we have 
replaced 5y~ with y~ and translations made to simplify 
the resulting expression). 



d'^Wtt" 



dH^ 



P0....,P„W^0 



1 r 



n! 



2N,'^ J (2^)3 

^e^P--v-{p\F-+-i^y-^y^)F\UQ)\p)- 



(45) 



In the next section, a re-summation of all such contribu- 
tions of arbitrary order n will be carried out to calculate 
the leading differential distribution of a hard jet in trans- 
verse momentum as a function of the length traversed in 
the medium. 
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V. RE-SUMMATION AND TRANSVERSE 
MOMENTUM DIFFUSION 

In the preceding section, the leading length enhanced 
twist-2n contribution was evaluated. As is obvious from 
Eq. (|45p . the essential combination of factors is that in- 
cluded in the square brackets. When this factor is not 
very small, all such terms (for all n) need to be re- 
summed into the generalized differential hadronic tensor, 
which will include contributions at all- twist. 

To simplify the re-summation, one further approxima- 
tion is required. In the preceding section, a model of 
the nucleus as a weakly interacting homogeneous gas of 
nucleons was used. The formal expression of this as- 
sumption is hidden within the dimensionful parameter 

a' 



The precise evaluation of such combinatoric 
coefficients is rather complicated, even for the case of 
next-to-leading twist [lO, [2l| . From general dimensional 
arguments, in the case of non-interacting nucleons, this 
factor may be approximated as. 



In the above equation, xt = {p\)/{2p ' q 



o\)/{2p+q ), where {p\) 
is the average transverse momentum that a nucleon may 
impart to a hard parton passing through it and G{xt) 
is the longitudinal gluon distribution function of the nu- 
cleon at XT, where 



G{x) 



dy- 



-ixp ' y 



2tt 



xp^ 



.{p\F^+'^{y-)F\+mp). (50) 



The scale of the gluon distribution is essentially the same 
as that used for the scatterings as well as for the final 
fragmentation i.e., m^. 

Following the form of Eq. the re-summed differ- 
ential hadronic tensor may be expressed in terms of a 
product of the leading order, leading twist, integrated 
hadronic tensor and a length dependent transverse mo- 
mentum distribution function (f>(L~ ,l±), 



~d?r~ 



(51) 



ci 



p 

2p+ 



(46) 



where, p is the nucleon density inside the nucleus and 
l/2p'^ originates in the normalization of the nucleon 
state. The remaining unknown coefficient is now 
considered to be independent of the order n. Straightfor- 
ward substitution the above approximation into Eq. (j45p . 
leads to a simplified form for the 2ii}^ order contribution 
to the differential hadronic tensor. The full hadronic ten- 
sor, which includes contributions from all orders, may be 
expressed as. 



i=0 



dH 



(47) 



While our analysis has focused on the hadronic tensor, 
it should be pointed out that W^'^ may be immediately 
supplanted with the differential cross section to produce 
a jet with a transverse momentum using Eq. In 
a sense, the entire process of the production and sub- 
sequent propagation of the hard quark through the nu- 
cleus has been factorized into two parts: C^W^" rep- 
resents the initial production in a hard scattering with 
a virtual photon and (j}{L~ ,l±) represents its propaga- 
tion and transverse momentum broadening. The entire 
length {L~) and transverse momentum (Zj^) dependence 
is contained entirely within the factor [(f){L~ Dif- 
ferentiating Eq. (|48|) with respect to L~ suggests the ob- 
vious relation for the differential transverse momentum 
distribution, 



Note that while i runs over all integers, each term refers 
to an even number (2i) of gluon insertions. Using the sim- 
plifications mentioned above and Eq. for the 2n"^ or- 
der term, we obtain the generalized differential hadronic 
tensor, 



(48) 



where, d'^W^'^ /d^l±^ is given by Eq. ([7]), and the constant 
D is given as. 



D 



d^yd'^p±_ 
2Nc J (27r)32p+ 



{p\F-^'^{y)F\^mp). 



X exp <; -I \ i2^V -p^-yi- 
'-PXtG{xt)- 



AN, 



(49) 



d<l>{L-,l^) 



= DV^i^4>{L-,li_), 



(52) 



which is a two dimensional diffusion equation in the vec- 
tor with the trace of the diffusion tensor given by 
Eq. (|49p . The initial condition is discerned from Eq. (O 



(53) 



The two dimensional diffusion equation, with the above 
mentioned initial condition, has the well known solu- 
tion [H, 



■ exp 



ADL- 



(54) 



Using the above diffusion equation, all moments of the 
transverse momentum distribution may be calculated. 
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The first non-zero moment is the average squared trans- 
verse momentum distribution at a length L~ , 



{IDl- = j d'li_ll4>{L- ,1^) = ADL- . (55) 

As a result, the transport coefficient of energy loss defined 
as the average transverse momentum squared gained by 
the jet per unit length traversed in the medium is given 
as (note: L~ = 2L the actual length traversed), 



similar result as the rigorous quantum mechanical anal- 
ysis performed above. In this case, one computes the 
propagation of a colored particle as it propagates under 
the action of fluctuating color fields which have a short 
color correlation length. This last condition is meant to 
mimic the effect of the color confinement at nucleon like 
distances as used in the preceding derivation. 

Imagine, as in the preceding section, that a colored par- 
ticle moves under the influence of a color Lorentz force in 
the (— ) direction. The equation of motion for a classical 
particle is, 



Nr. 



pxTG{xT,m'^)\^^o, (56) 



It should be pointed out that in an earlier work where a 
diffusion equation was used to understand the systemat- 
ics of transverse momentum broadening of jets (23| . the 
diffusion tensor D was defined to be four times that in 
Eq. (|49|). such that one obtains, the mean square trans- 
verse momentum broadening as (/^) = DL^ and q = 2D. 

The final numerical value of the scale of the gluon dis- 
tribution function deserves some discussion. We have 
computed a class of all twist corrections to what amounts 
to a leading order (LO) hard scattering process. In the 
numerical implementation of such processes, there exists 
only one large scale i.e., Q^. Thus, while in principle 
the multiple scattering and final fragmentation occur at 
a softer scale, in LO calculations one may set = 
for both the gluon distribution as well as the final frag- 
mentation. If the hard scattering were computed at NLO 
this may no longer be the case. 



(57) 



where, t is the time in the observers rest frame, Q'^{y) is 
the color charge expectation of the particle at location y, 
F"-^^ {y) is the color field strength at that location and 
V is the three velocity written in four component form 
i.e., V = [Ij'S] in time-space components. In light cone 
components the vector v is given as, v =[[1 + Uz)/2, 1 — 
v±\ . Note that for a light like particle + v'j_ = 1, 
thus V is akin to a unit vector. For a particle moving 
predominantly in the negative z direction as in the case of 
the preceding sections, Vz < and \vz\ » \v_l\. Hence, 
one may make the approximation, 



F^.v-' ~ 2F+. 



(58) 



As a result, the squared transverse broadening in mo- 
mentum experienced by the light-like particle is given as 



VI. LORENTZ-LANGEVIN ANALYSIS AND 
CLASSICAL PROPAGATION 

The field theory calculation presented in the preced- 
ing sections demonstrates that, in the limit of no radia- 
tion, the transverse dynamics of hard collincar jets pass- 
ing through dense matter may be understood in terms 
of a transverse momentum diffusion equation. As such, 
no clear physical picture has emerged from all the ap- 
proximations carried out. The situation at this point is 
not dissimilar to that encountered in the hard thermal 
loop (HTL) effective theory of QCD [H, (H] where hard 
momenta {p ~ T) were integrated out to find effective 
propagators for soft modes (p ~ gT). Identical results 
were also derived if one assumed the more physical pic- 
ture of the hard modes to be classical particles moving 
under the influence of a soft color field which represented 
the soft modes [l^ . In this section, we demonstrate that 
an almost identical situation is true for the case of higher 
twist jet broadening. 

In the following, an alternative calculation of the trans- 
verse broadening of an energetic jet as it passes through 
a dense colored medium will be carried out classically. It 
is not a complete surprise that such a derivation yields a 



\^p±iy)\ 



(59) 



y 

= g'J dy^dy^Q^Q''F\+{y^)F''^"{y^), 

3/0 

where, the substitution y~ = 2t has been made for a light 
like particle traveling in the negative z-direction. Taking 
the expectation of the product of gluon field strength 
operators, one may now institute the condition that the 
medium has a short color correlation length, i.e.. 



dy^dy^{F\+{y^)F''^ (y^)) 
J dy^{F\+F''-'')y-^^f. (60) 



The average over the color of the gluon field correlators 
may now be carried out (ignoring position coordinates). 



{F'^F'') 



nab 



{N^ - 1) 



{F^F" 



(61) 



The color charges, Q"" represent the expectation of the 
color matrices in the state of the colored parton. As a 
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result 5°''^Q°'Q^ = Cp for a quark parton. Hence, one 
obtains the mean square transverse momentum after a 
parton travels a distance from i/q to y~ as given by, 



Ap_L(y) 



27r 



(62) 



£_ j pa+a 

2 



-y_ 

2 



Using the same normalization as in the preceding sec- 
tion for turning a nuclear state into a gas of nucleon states 



I.e. 







2p+ 



{p\0\p). 



one obtains identical expressions for the mean square 
transverse momentum per unit length as obtained from 
the higher twist analysis. This leads to the contention 
that the re-summation of length enhanced twist correc- 
tions on to an outgoing parton line is effectively described 
as the classical propagation of hard colored partons in 
soft fields, as long as the correlation lengths of the soft 
fields are short. 



VII. DISCUSSIONS AND CONCLUSIONS 

The propagation and energy loss of hard partons in 
dense matter is a topic of much current interest. In the 
current effort, a detailed analysis of the transverse mo- 
mentum diffusion of jets produced in DIS off large nuclei, 
as they pass through cold nuclear matter was carried 
out. Radiation and radiative energy loss were ignored, 
though, in reality, such partons will tend to radiate and 
lose a substantial fraction of their energy as they traverse 
the dense matter. Complications arising from energy loss 
were ignored in an effort to isolate the systematics of the 
propagation of a single parton in a fluctuating color field 
without radiation. Radiation will be included in a sub- 
sequent effort to complete the formulation of the propa- 
gation of hard partons through extended dense matter. 

The leading length enhanced higher twist corrections 
which produce maximal broadening were identified and 
re-summed onto the out going parton line. The final 
transverse momentum distribution at a given length tra- 
versed was shown to obey a two-dimensional diffusion 
equation. The diffusion equation was solved and the 
mean square transverse momentum per unit length (g) 
calculated. This is the first effort, to the knowledge of 
the authors, which presents a first principles derivation of 
the transverse momentum distribution function (within 
the higher twist formalism) and a relation between the 
diffusion tensor and the transport coefficient Our re- 



sults differ from those of an earlier attempt [ll| to re- 
sum length enhanced multiple scattering diagrams, which 
found a shift, rather than a broadening, of the transverse 



momentum distribution. We believe that our result is 
more broadly applicable. 

The results for q obtained from the higher twist analy- 
sis were shown to be identical to that obtained for a point 
like colored particle moving in a fluctuating color field if 
the color correlation length of the fields is assumed to be 
short. This places the effective picture of higher twist re- 
summation on the same footing as that of the HTL effec- 
tive theory which also admits an almost similar picture. 
This suggests the possibility of constructing an effective 
field theory applicable to the propagation of hard jets in 
soft dense matter. 

The derivations in this effort have been restricted to 
parton propagation. The construction of an effective the- 
ory will require also the use of re-summed vertices. These 
will be derived in a subsequent publication when radia- 
tive processes will be dealt with. The inclusion of such 
contributions will allow for an all-twist calculation of ra- 
diative energy loss and comparisons with the energy loss 
formalism of Ref. The inclusion of radiative energy 
loss in an all twist formalism will allow, for the first time, 
a means to study the evolution of virtuality of the hard 
partons as they originate in hard collisions, traverse dense 
matter, escape and fragment into hadrons. 

Confrontation with experiment f27| will require a con- 
volution of the above formalism with both single [J and 
multi-hadron [28l . [29| fragmentation functions. Even 
without the inclusion of energy loss, the above formal- 
ism may be directly applied to the measured transverse 
broadening of hard jets in cold nuclear matter as mea- 
sured by the HERMES experiment at DESY Such 
comparisons will be carried out in an upcoming publica- 
tion. 



VIII. ACKNOWLEDGMENTS 

The authors wish to thank R. J. Fries and X.-N. Wang 
for helpful discussions. This work was supported in part 
by the U.S. Department of Energy under grant DE-FG02- 
05ER41367. 



APPENDIX A: GLUON FUSION DIAGRAMS 

In this appendix, we briefly sketch the calculation of 
diagrams that include gluon fusion contributions such as 
in the right panel of Fig. |4l Such diagrams originate 
from the inclusion of interaction terms from the QCD 
Lagrangian density, of the form. 



^ (d^A^y,) - d^A^y,)) r''^A'^'^{y,)A''^{y,). (Al) 

A term in the expansion of the hadronic tensor that 
would contain a fusion contribution will have the gen- 
eral form [ignoring (— ) signs and factors of i and g], 
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M^aT = J dydy,dyg{J^{y)4>{ypWt''A'i^{yp)4>{yp) (A2) 

x(5„A^(y,) - d^A^y,)) r''A'^^{y,)A'\y,).r). 

In the above equation. and J" are the currents at 
the locations y and which couple to the photon and 
there is a contraction between the gluon vector potential 
appearing in the first line with any of the three gluon vec- 
tor potentials in the second line, this contraction repre- 
sents the one gluon propagator in the left panel of Fig. H) 
There are multiple reasons why such contributions are 
small. One reason is that unlike the ladder diagrams of 
Fig. O the introduction of an extra coupling constant is 
not accompanied with the introduction of an extra length 
integral. As the momentum in the one gluon propagator 
is set by requiring that the two quark lines it connects to 
go nearly on shell, the gluon propagator cannot itself go 
on shell and as a result there is no length enhancement 
from the gluon propagator. 

This contribution is also small due to color confine- 
ment and the condition that the forward momenta of all 
partons originating in the nucleus are not small enough 
to be in the saturation regime. The contraction of one 
of the vector potentials located at yg [in the second line 
of Eq. (|A2[) ] with that at j/p, requires that the remain- 
ing two vector potentials be contracted into the nuclear 
state. As both vector potentials are at the same location, 
they must be contracted into the same nucleon. How- 
ever, due to Eq. (|1T|). this requires both gluons to be in 
a color singlet state. The presence of the antisymmetric 
factor /"'"^ renders such contributions vanishing. Dia- 
grams with three gluon vertices produce non-vanishing 
contributions if another gluon (at location say y') were 
to be contracted with the two gluons (which fuse into 
one gluon), into the same nucleon state. Such contribu- 
tions automatically require that another space integral 
over y' be restricted to be in close proximity to yg. Also, 
such contributions necessarily require the introduction of 
the expectation of 3-parton operators in a nucleon state, 
which, as we have argued previously, are small compared 
to the expectation of 2-parton operators outside the sat- 
uration regime. 



APPENDIX B: NON-CENTRAL CUT DIAGRAMS 

In Sect. IV, the derivation of the transverse momen- 
tum diffusion of a hard parton was isolated to the case 
of a central cut. In Eq. pQ]) and the ensuing discussion, 
this was identified as the dominant contribution to the 
differential hadronic tensor at a given order. In this ap- 
pendix, we focus on diagrams at order n + n' = 2m in 
terms of the strong coupling constant g, where n ^ n' . 
In particular, we focus on the case where n = m — 1 
and n' = m + 1. The starting point of this analysis is 
Eq. (PT|) . As in Sect. IV, the hard part is expanded in 



a Taylor series in transverse momentum as in Eq. (|36p . 
resulting in the integrand (written without phase factors 
and light-cone time ordering 0-functions as), 



H{p^,q ,p±,y)S{pj_,y) 

m—1 m+1 ^ 

n I[izr^Hip+,q-,p^,p'i) 



t\ f=\ 9p\ dp\ 



(Bl) 



As in the case of the central cut, the 0-functions 
time order the gluon insertions in incremental order 
from the photon vertex to the cut line. Each occur- 
rence of the factor e^^ ^'^ ^p\^A^ [y"^) may be converted, 
using integration by parts, to —ie'^P^'^^ d]_^A'^ {y"^) ~ 
-F±^{y^). As a result we obtain an over all sign 



-le 



of (— i)*^™ — —1 compared to the central cut. 




FIG. 6: An order 2m, contribution to VF*"^ from a non-central 
cut 

The nuclear state is decomposed as before into a 
pseudo-free gas of nucleons and pairs of gluon insertions 
in singlet combinations is contracted into the m chosen 
nucleons. Unlike the central cut, there is an excess of 
two gluon insertions on the left hand side of the cut as 
shown in Fig. [6l where the rectangular blobs represent 
single nucleon states. This extra pair is contracted into 
the same nucleon. For this specific example, the m*'' and 
{m + 1)*'^ gluon on the left hand side will be contracted 
into the same nucleon, as required for maximum length 
enhancement Contractions within the same nucleon, due 
to color confinement, tend to constrain the locations of 
the field operators F^''{y^) and F^^(y„+i) to a region 
y7n+i ~ Dm 1^ Vconf the Confining distance. Thus con- 
tractions between non-adjacent gluons, due to light-cone 
ordering 0-functions will constrain the length integrations 
of all the intermediate gluons to be within the size of a 
nucleon. 
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The transverse momentum integrations may now be 
carried out with the observation that in a very large nu- 
cleus, two point correlators only depend on the relative 
transverse distance between the two insertions. The final 
result is similar to that in Eq. ((34)) for i < m — 1 [see 
Eq. p4|) ]. The resulting (5-functions for the case of the 
m"^ and (m + 1)"^ gluon which are both on the left of 
the cut, are slightly different, i.e., 



dyTdSy 



(m+l) 



2 e2 (rn , (r?i + l) 



dSy'le 



X (t'^-Fl-+J6y^/2)t^-F''^-^'"\-5y'l/2)). (B2) 

The transverse momentum ^-function in the above equa- 
tion imposes the condition that the transverse momen- 
tum brought in by one gluon is taken out by the other. 
Thus a pair of adjacent gluon insertions, in the limit of 
transverse position invariance, infuses no net transverse 
momentum, into the final cut line. 

Incorporating the m transverse momentum delta func- 
tions, leads to the following simplified expression for the 
hard part of the integrand. 



2 




(B3) 



where, (5y„, = y„, -y^^^+iy 

One notes the absence of the two transverse momenta 
pT'Pj™"''^^ from the argument of the i5-function which 
determines the momenta of the cut- line This will 

also be the case if the extra gluons were instead consid- 
ered on the right hand side of the cut line. The essential 
point is that in the case of translational invariance in a 
large nucleus, the leading length enhanced contribution 
from double scattering off a single nucleon produces 
no net transverse momentum broadening. As a result, 
such diagrams contribute to terms with a lower number 
of transverse momentum derivatives. The transverse 
momentum derivatives that act on the phase factors lead 
to suppressed contributions. An evaluation of the double 
derivatives on the phase factors, in the limit of vanishing 
transverse momenta, will leave an over all factor of 
iSy^/{q~). As Sy^^ is restricted to the size of the 
nucleon and the entire contribution is suppressed by a 
factor of q~ , such terms arc paramctrically suppressed in 
the evaluation of the transverse momentum distribution. 
Such terms are however important in terms of unitary 
corrections to the total integrated cross section as will 
be demonstrated in the next appendix. 



APPENDIX C: UNITARITY CORRECTIONS 

In the evaluation of the transverse momentum distri- 
bution [(j){L~ , l±)] in Eq. as given by the solution of 
the diffusion equation [Eq. ([52]) ]. the focus was limited to 
the central cut diagrams as in Fig. [5l It was argued that 
diagrams with non-central cuts do not contribute to the 
leading behavior of the transverse momentum distribu- 
tion. These diagrams, however, play an important role 
in the unitarization of the integrated cross-section. The 
unitarity requirement states that 



dJd'^l±(j)(L-,lj_) 



= 0. 



(CI) 



One notes that the above requirement is satisfied by 
Eq. (j54p . The unitarity requirement was tacitly assumed 
in formulating the overall normalization factor of the 
Gaussian distribution in Eq. In what follows, we 

discuss the origin of such a constraint. 

The starting point is that of Eq. pT|) . In the interest of 
simplicity we focus on the case oin+n' = 2. Cases with a 
higher number of gluon insertions follow a similar, albeit 
a more complicated pattern. For the case of two gluon 
insertion there arc three cuts: the central cut with n = 
n' = 1, the left cut with n — Q,n' = 2 and the right cut 
with n = 2, n' = 0. The contribution to the central cut 
may be immediately evaluated from Eq. (|3ip . by setting 
n ^ n' = I. Integrating out the transverse momentum of 
the cut-line 1 1 we obtain the hadronic tensor as 



= f f[dy-dY^f[dy'^d 

^ ^ — n A — 1 



- 72 n 

y L 



i=0 3 = 1 



^ ^2„i ^2^/0 



(27r)2 (27r)2 



ix]jp+y. ip\-y\^ 



1=0 



X He^^'^P^y'J e-^P'Wie{y^ -y-)0iy';) 

X {A;p\iiy-,y^)j+^{Q)Tr [eA+{y^,yi) 
X (C2) 

The nuclear state is now decomposed into a gas of nu- 
cleon states with the quark wave function insertions re- 
stricted to one nucleon and the the gluon vector potential 
insertions restricted to another nucleon. 

In what follows we focus solely of the matrix element 
of the two gluon insertions in a nucleon, which is referred 
to as the soft part S, 
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= gi{p±-p'±)-y±gi{p±+p'±)-Sy±/2 

X {p\Ai{6yj2)A+i-5yj2)\p), (C3) 

where, the transverse positions have been recast as y + 
Sy/2 and y — 5y/2. Transverse position invariance in 
a large nucleus is invoked in the last line of the above 
equation where the two point correlator is set to be inde- 
pendent of the mean transverse location y±. The integra- 
tion over the mean transverse location can be carried out 
to obtain the two dimensional delta function (5-function 
5'^{p]_ — The hard part, which essentially consists 

of phase factors may now be expanded in a Taylor se- 
ries in the transverse momenta p±,p'\ as in Eq. ((36)) . 
Carrying out the derivatives and setting the remaining 
transverse momenta to zero, gives the part that depends 
on the gluon insertions as 

G ^ I dy^dy'-e{y^ -y-)e{y'-) 

X ^ ^^\^^'\ p\F^{yi)F+"iy'V)\p)- (C4) 

In the above equation, due to color confinement, y^ — y'l 
is always restricted to be smaller than the size of a nu- 
cleon. However, the mean location +y'i )/2 is unre- 
stricted and as a result this length integration provides 



an overall factor of and is the origin of the length 
enhancement of the part of the cross section which cor- 
responds to the central cut. 

Setting n = 0, n' = 2 and n = 2,n' — 0, and following 
an almost similar method as above for the central cut 
yields the contributions from the left and right cuts. As 
shown in the preceding appendix, the terms originating 
from the left and right cuts provide almost identical con- 
tributions, except for an over all sign and the argument 
of the ^-functions. Relabeling the locations: e.g., y2 = y'l 
in the case of the left cut (and similarly for the right cut 
to have the same set of position labels for all three cuts), 
yields the combination of ^-functions as 



e = 0{y- -y-)e{y\)-eiy\ -y^)e{y^-y-) 
- e{y^ -y'-)0{y'-). (C5) 

As has been demonstrated earlier 0, [l^, [ll|, [l^, the 
above combination, along with the constraint from con- 
finement destroys the length enhancement by restricting 
the average longitudinal location (yf +y'i )/2 to also lie 
within the size of a nucleon. The combination of all three 
cuts, in the evaluation of the integrated distribution, thus 
removes the overall factor of L~ and as a result the inte- 
grated distribution obeys the condition of Eq. (jCip . The 
general case with 7i-gluon insertions may also be demon- 
strated to contain such unitarity corrections by a similar 
method of combining contributions from all cuts. 



[1] K. Adcox et al. [PH ENIX Collaborat ion], Nucl. Phys. 
A 757, 184 (2005) arXiv:nucl-ex/0410003 ; J. Adams et 
al. [STAR Collaboration], Nucl. Phys. A 757, 102 (2005) 



arXiv:nucl-ex/0501009 



[2] M. Gyulassy and M. Plumer, Phys. Lett. B 243, 432 
(1990); X. N. Wang and M. Gyulassy, Phys. Rev. 
Lett. 68, 1480 (1992); R. Baier, Y. L. Dokshitzer, 
A. H. Mueller, S. Peigne and D. Schiff, Nucl. Phys. B 483, 
291 (1997) arXiv:hep-ph/9607355 ; Nucl. Phys. B 484, 
265 (1997) arXiv:hep-ph/9608322 ; B. G. Zakharov, 
JETP Lett. 63, 952 (1996) arXiv:hep-ph/9607440]; 
JETP Lett. 65, 615 (1997) arXiv:hep-ph/9704255'; 
Phys. Atom. Nucl. 61, 838 (1998) [Yad. Fiz. 61, 924 
(1998)] arXiv:hep-ph/9807540 . 

P. Arn old, G. D. Mo ore and L. G. Yaffe, JHEP 0206, 030 
(2002) (arXiv:hep-ph /0204343 ; S. Jcon and G. D. Moore, 
Phys. Rev. C 71, 034901 (2005) arXiv :hep-ph /0309332 ; 
S. Turbide, C. Gale, S. Jeon a nd G. D. M oore, Phys. Rev. 
C 72, 014906 (2005) tarXiv:hep-ph /050 2248 . 
[4] X. F. Guo and X. N. Wa ng, Phys. Rev. Lett. 
85, 3591 (2000) arXiv:hep-ph/6005044 ; X. N. Wang 
and X. F. Guo, Nu cl. Phys. A 696, 788 (2001) 
arXiv:hep-ph/0102230| ; B. W. Zhang and X.-N. Wang, 



Nucl. Phys. A 720, 429 (2003) arXiv:hep-ph/0301195^; 
[5] M. Gyulassy, P. Levai and L Vitev, Phys. Rev. Lett. 85, 
5535 (2000) [arXiv:nucl-th/0005032| ; Nucl. Phys. B 594, 



[6] 
[7] 
[8] 
[9] 

[10] 

[11] 

[12] 

[13] 
[14] 
[15] 

[16] 



371 (2001) [arXiv:nu cl-th/0006010| ; 
U. A. Wiedemann 



Phys." B 588, 303 (2000) 



arXiv: hep-ph/0 302184, . 

Gyulassy, Nucl. Phys. A 733, 265 



Nucl. 

arXiv:hep-p h/00051 29 . 

C. A. Salgado ahd U. A. Wiedemann, Phys. Rev. D 68, 

014008 (2003) 

M. Djordjevic and M. 
(2004) arXiv:nucl-th/0310076 . 

V. N. Gribov and L. N. Lipatov, Yad. Fiz. 15, 781 (1972) 
[Sov. J. Nucl. Phys. 15, 438 (1972)]; Yu. L. Dokshitzer, 
Sov. Phys. JETP 46, 641 (1977) [Zh. Eksp. Teor. Fiz. 
73, 1216 (1977)]; G. AhareUi and G. Parisi, Nucl. Phys. 
B 126, 298 (1977). 
M. Luo, J. w. Qiu and G 
377 (1992); M. Luo, J. w 
Rev. D 50, 1951 (1994). 
R. J. Fries, Phys. Rev 
X. 



Sterman, 
Qiu and 



Phys. Lett. B 279, 
G. Sterman, Phys. 



D 68, 074013 (2003) 
Guo, Phys. Rev. D 58, 



'arXiv:hep-ph/02 09275, 

114033 (1998) ar Xiv:hep-ph/9804234 
J. Casaldcrrcy-Solana and D. Teaney, JHEP 0704, 039 
(2007). arXiv:hcp-th/0701123 . 

J. w. Qiu and G. Sterman, Nucl. Phys. B 353, 105 (1991). 
J. w. Qiu and G. Sterman, Nucl. Phys. B 353, 137 (1991). 
J. C. Collins, D. E. Soper, and G. Sterman, in Pertur- 
bative Quantum Chromodynamics edited by A. Mueller, 
(World Scientific 1989), and references therein. 
J. Casalderrey-Solana and X. N. Wang. arXiv:0705.1352l 



17 



[hep-ph]. 

[17] A. H. Mueller, Nucl. Phys. B 335, 115 (1990); Nucl. 

Phys. B 415, 373 (1994). 
[18] L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 

2233 (1994) [arXiv:h ep-ph/9309289 ; Phys. Rev. D 49, 

3352 (1994) [arXiv:h ep-ph/9311205 . 
[19] X. Ji, J. W. Qiu, W. Vogelsang and F. Yuan, Phys. Lett. 

B 638, 178 (2006) arXiv: hep-ph/0604 128 . 
[20] A. Majumder, E. Wang and'~X7'NrW ang, Phys. Rev. 

Lett. 99, 152301 (2007) [arXiv:nucl-th/041206i; . 
[21] J. Osborne and X. N. Wang, N ucl. Phys. A 710, 

281 (2002) 'arXiv:hep-p h /0204046| ; A. Majumder and 

X. N. Wang, to appear. 
[22] L Sneddon, Elements of partial differential equations, 

McGraw HiU, Singapore (1957). 
[23] A. Majumder, B. M uller and S. A. Bass, Phys. Rev. Lett. 

99, 042301 (2007) tarXiv:hep-ph/ 061 1 135, . 
[24] E. Braaten and R. D. Pisarski, Nucl. Phys. B 337, 569 

(1990). 

[25] J. Frenkel and J. C. Taylor, Nucl. Phys. B 334, 199 



(1990). 

[26] J. P. Blaizot and E. lancu, Phys. Rev. Lett. 70, 3376 

(1993) arXiv:he p-ph/9301236 ; Nucl. Phys. B 417, 608 

(1994) (arXiv:hep -ph/9306294 . 

[27] K. Adcox et al. [PHENIX Collaboration], Phys . 

Rev. Lett. 88, 022301 (2002) [arXiv:nu cl-ex/0109003 | ; 

C. Adler et al. [STAR CoUaboraSTo^JT^Phys. Rev. Lett. 

89, 202301 (2002) arXiv: nucl-ex/0206011 
[28] A. Majumder an d X7N7Wa ng, Phys. Rev. D 70, 014007 

(2004) arXiv:hep-ph /0402245t. A. Majumder, J. Phys. G 
30, S1305 (2004) ar Xiv:hep-ph /0404292 . A. Majumder 
and X. N. Wang, jaxXivdiep-ph/0410078 

[29] A. Majumder and X. N. Wang, Phys. Rev. D 72, 034007 

(2005) arXiv:hep-ph76411174 . 

[30] A. Airapetian et al. [HERME S Collabora tion], Eur. 
Phys. J. C 20, 479 (2001) [arXiv:hep-ex /0012049l; 
A. Airapetian et al. [ HERMES CoUaboratio n] , Phys. 
Lett. B 577, 37 (2003) [arXiv:hep-ex/0307023|. 



